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ABSTRACT 

Maxwell's equations and the equations governing charged particle dynamics 
are presented for a rotating coordinate system with the global time coordinate of 
an observer on the rotational axis. Special care is taken in defining the relevant 
entities in these equations. Ambiguities in the definitions of the electromagnetic 
fields are pointed out, and in fact are shown to be essential in such a system of 
coordinates. The Lorentz force is found to have an extra term in this frame, 
which has its origins in relativistic mass. A related term in the energy equation, 
which allows inertia to be gained even during strict corotation, suggests ways 
existing pulsar magnetosphere models may be modified to match observed 
'braking indices' more closely. 

Subject headings: magnetic fields — relativity — pulsars: general 
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1. Introduction 

In pulsar magnetosphere theory and other applications (Chedia et al. 1996; Cheng 
1984; Fawley, Arons, & Scharlemann 1977; Hones & Bergeson 1965; Schiff 1939), it is 
often convenient to adopt a rotating coordinate system as a frame of reference. Such a 
frame is non-inertial, and has the usual Coriolis and centrifugal fictitious forces present as a 
well-known consequence. However, electrodynamic processes are often of primary interest, 
and this leads to several difficulties. In particular, the interpretation of various effects and 
even the definitions of the electric and magnetic fields are necessarily ambiguous (some 
consequences of this were noticed by Backus 1956). It is the purpose of this work to explore 
these matters in some detail and thereby provide clarification. In particular, we will first 
present Maxwell's equations and then the Lorentz force in a system of rotating Cartesian 
spatial coordinates, with the global time coordinate of a (stationary) observer on the axis 
of rotation. 

As will be seen, it is difficult to define the electric and magnetic fields in this frame. 
This situation worsens with distance from the axis of rotation, and becomes critical at the 
'light cylinder' distance, tl = 1/^, where Q is the angular frequency of rotation (here and 
throughout, we use units where c = 1). This situation will be shown to be an essential 
feature of all such frames whose metric tensor g^ u has off-diagonal elements. For this reason, 
all electromagnetic quantities used here will be given careful definitions, with reference to 
their values in an inertial non-rotating frame (where the ambiguities disappear). 

Among the most surprising consequences which will be shown are that the Lorentz 
force acquires an extra 'relativistic mass' term in this frame, and particles may gain inertia 
even when undergoing strict corotation. Although negligible at low altitudes, these effects 
also can become critical at altitudes approaching the light cylinder. In pulsar models, it is 
precisely in this regime where the magnetic polar magnetospheric currents are expected to 
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close; the details of this process remain an open question whose answer is of key importance 
for the models. It is therefore hoped that this paper will serve both as a reference for those 
who wish to adopt rotating coordinates when studying electrodynamic processes, and as a 
guide for interpreting the various phenomena in these coordinates. 

2. Definitions 

We must first define our coordinate systems; the remainder of the mathematics is then 
straightforward tensor algebra. To avoid an overabundance of 'primes,' all coordinate- 
dependent quantities with 'primes' attached shall refer to the inertial (non-rotating) frame, 
and 'unprimed' quantities will be used for the rotating frame. Greek letters will be used for 
indices which range over all four spacetime components, for example, U a = ([/*, U x , U y , U z ), 
while Latin letters will be used when variation is to be made over the three ordinary spatial 
indices (x,y,z). The usual summation convention is adopted. 

We start with an inertial nonrotating Cartesian coordinate system (£', x', y', z') with 
the 'flat' (ignoring gravity) Minkowski metric 

(ds') 2 = -{dt'f + (dx') 2 + (dy'f + (dz'f (1) 

where ds' is the infinitesimal proper distance for a spacelike interval (and is coordinate- 
independent, although the 'prime' is left on for clarity). If the interval is instead timelike, 
we define its proper time dr' (also coordinate-independent) by replacing (ds') 2 by — (dr') 2 
above. We write the above as 

(ds 1 ) 2 = g^vidx^' dx v \ (2) 

defining the metric tensor g^> u '. 

Next, we define coordinates which are rotating with respect to these, with angular 
velocity f2 counterclockwise about the z' axis. Thus, t = t',z = z',x = x' cos Vlt' + y' sin fit', 
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and y = y' cos fit' — x' sin fit'. Inverting gives x' = x cos fit — y sin fit, and 

y' = y cos fit + x sin fit. Note that the time coordinate is identical to that in the 

nonrotating frame: we are not measuring time as observers rotating with the frame would 

(which we shall see has its advantages as well as disadvantages). Naturally, we will have an 

'ergosphere' beyond r = r L , where all material particles must have velocities opposing the 

rotation. 

The differential coordinate transformation is dx^ = t^ v ,dx v \ which defines 
the coefficients A^, = dx^/dx u '. The inverse transformation defines the coefficients 
A£' = dx»'/dx u , so that A£, R"J = 6£. 

The metric g^ u > transforms as a (symmetric) second-rank covariant tensor: 

9llv = kfflgjp (3) 



which gives here 



/-(1-fiV) -toy fix 0\ 
-fly 1 

fix 10 



(4) 



V o 001/ 

(Although the use of cylindrical coordinates would obviously simplify some expressions, 
such as this one for g^ u , we stick to rotating Cartesian spatial coordinates to avoid the added 
confusion cylindrical coordinates give to covariant vs. contravariant components of tensors.) 
Notice that the coefficient of dt 2 vanishes on the cylinder fl 2 g 2 = 1, where g = \J x 2 + y 2 
is the distance from the rotation axis (not to be confused with the charge density p). 
However, the determinant of is always —1 due to the off-diagonal 'space-time' terms, so 
there is no mathematical difficulty in using these coordinates beyond g — tl- This would 
not be true if we had picked the times measured by local corotating observers as our time 
coordinate (as in Chedia et al. 1996), which would also have the undesirable property of 
having the definition of t vary with g. 



- 6- 



We shall need the inverse of g^ v . It is 
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The presence of nonzero off-diagonal terms in the metric tensor and its inverse remind us 
that these coordinates are not an orthogonal system, so that vectors and their associated 
one-forms (covectors) do not generally 'point' in the same direction. This will be shown 
to ultimately be the reason why the notion of separate 'electric' and 'magnetic' (spatial) 
vector fields runs into serious trouble in this frame. 



3. Maxwell's Equations 



The relativistically covariant way to discuss electrodynamics in an arbitrary coordinate 
system is to introduce the antisymmetric second-rank tensor F which has as covariant 
entries in the 'primed' (nonrotating) frame (Misner, Thorne, & Wheeler 1973; Weinberg 
1972; Wald 1984) 
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where E 1 ' are the (inertial-frame) components of the electric field E, B 1 ' are the components 
of the magnetic field B, and we have assumed vacuum permeability and permittivity. The 
components of F^ u in any other frame are then found in the same way as for g^ v in equation 
(3): 

(7) 
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We find 



F = 



( -E* 
E x 
E y -B z 



-E y —E z \ 
B z -B y 
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(8) 



\E Z B y -B x / 
where E l and B l are respectively defined as the projections of the electric and magnetic 

fields (as measured in the inertial frame) onto the new (contravariant) spatial coordinate 

vector basis. For example, B = B l 'e v = B l ei (at corresponding locations of the two 

coordinate systems; we are abusing notation), where the are the (spatial) contravariant 

coordinate basis for vectors in the unprimed coordinates. That is, = (d/dx l ), holding x J 

fixed for all j ^ i. We shall treat the E 1 and B % as numbers (rather than the contravariant 

components of vectors) in what follows, for clarity. F is the fundamental tensor of interest, 

while E and B are not parts of 4- vectors and so do not transform as such. 



We have also defined the vector 

E = E+(C2xr)xB 



(9) 



which, from equation (8), seems to be the electric field in our new coordinates. However, 
we will also need the contravariant formulation of F, easily found by raising the indices on 
equation (8) with the metric: 



where we have defined 



It is straightforward (and reassuring) to verify that the invariants 
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(11) 



B B - E E = -F^F^ 
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and 



E B = X -F» v *F flv 



are reproduced by equations (8-11), where *F IXV = l/2e fll/a/3 F a ^ is the 'dual' tensor of F, 
and is the totally antisymmetric fourth-rank tensor (4-form) with £0123 = 1- 

Notice that equations (9) and (11) do not represent a local Lorentz transformation 
to our rotating coordinates, as our frame is not locally Lorentz (due to the global time 
coordinate used). Clearly, the interchanges between the 'normal' and 'tilde' versions of E 
and B when passing from the covariant to the contravariant representations of F show a 
problem when trying to define the fields in these coordinates. We will see from Maxwell's 
equations that we have no obvious way of choosing E, E, or some combination of these 
objects as the 'electric field' in this coordinate system (for reasons to be detailed in section 
5). Some authors (Cheng 1984; Fawley et al. 1977; Hones & Bergeson 1965) simply follow 
the precedent set by Schiff (1939), who was resolving an interesting paradox involving 
rotating frames and Mach's principle, and who arbitrarily used the covariant choices E and 
B. Schiff correctly pointed out that it is the components of the mixed tensor F^ v which 
dictate particle motion; we shall return to this point after deriving Maxwell's equations for 
the rotating frame. 

Maxwell's equations in relativistically covariant form can be written (Weinberg 1972) 

S (V=^) =W ^ (12a) 



dx a 
and 

where g is the determinant of (simply —1 here), and J is the current density 4- vector. 
In the inertial frame, we have J' a = (p',j') which transforms to J a = (p, j), where p = p', 
and j = j' — p(f2 x r). Note that the convection of charge due to the rotation alters the 
effective current density. 
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Equation (12a) produces 



and 



while equation (12b) gives 



and 



V • E = 4vrp (13a) 



r?F 

VxB = 47rj + — , (136) 



V • B = (13c) 



VxE+ — = 0. (13d) 
at 



Note the appearance of both 'normal' and 'tilde' components of E and B. Equations (13) 
constitute Maxwell's equations in this frame. Here, V is the usual 3-dimensional gradient 
operator, if rotating Cartesian coordinates are used. 

We combine equations (13a) and (13b) to get an equation for charge conservation: 

| + V.J = (14) 

Equations (13) and (14) suggest that there is no difficulty in interpreting p and j as charge 
and current densities in this frame. However, J a is (unlike E and B) a genuine 4-vector, 
with covariant components J a = (— p(l — VL 2 q 2 ) + (ft x r) ■ j, j'). 



4. Geodesies and Lorentz Force 

To explore particle dynamics, we use the geodesic equation, modified by the Lorentz 
force when necessary. First, we find the Christoffel symbols from the metric: 

f m _ 1 n ^( d 9ua dg v p _ dg af3 \ 

L ^~2 9 [dx? + dx« dx»y 1 ' 

The only nonvanishing ones are found to be Y^ t = —Q 2 x, Y\ t = —fl 2 y, T* t = = —Q, and 



1 xt — L tx — 
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We now define the 4-velocity of a massive particle to be 

U" = ^- (16) 
dr 

where r is the proper time of the particle. Notice that this gives U^U^ = g fiU U fl U u — — 1, 
by the definition of proper time. A particle subject to no external forces travels along 
geodesies, which are the solutions of 

dU 11 TTVTjX 



, + T^ X U U U A = 0. (17) 
dr 

The temporal equation is trivial: 

^ = 

dr 2 ' 

giving 

dt 

— = constant = 7. 
dr 

We note that 

d dt d d 
dr dr dt ^ dt 1 

so that = (7,7V), where v = dx/dt. Plugging this into g ilv U ii \J v = —1 gives, after 
rearranging, 

7= 7rh/> (18) 

where v' = v + O x r, the particle's velocity as measured in the nonrotating frame of 
reference. Thus, 7 (which is conserved in the absence of external forces) is just the particle 
energy in the nonrotating frame (in units of mc 2 ). Note, however, that 7 is not equal to 



l/\/l — v 2 , which is fortunate since v > 1 is possible in this frame. 

The spatial geodesic equations give (after cancelling a common 7) 

dw 

- = -J]x(fixr)-2(Ox v). 

(XL 



- 11 - 



This is just the familiar expression for the Coriolis and centrifugal forces in a rotating 
frame. They are still valid relativistically, even beyond the light cylinder. 

When external forces are involved, we need only to introduce all forces (per unit mass) 
as 4- vectors f a , and insert them on the right hand side of the geodesic equations in a 
covariant manner: dU a /dr + T^ U U^U U = f a . 

The Lorentz force is given by the contraction of F^ u with the 4- velocity U a : 

f a = -9,uF a W (19a) 
m 

or equivalently by using the 'mixed' representation of F: 



dr m 



+ K X U"U' 

where q is the charge on the particle, and m its mass. The mixed tensor is here 
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Substituting this Lorentz force into the geodesic equations gives 

v + (ft x r) 



dr^ = q_ 
dt m 
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m 



for the temporal equation, and 
d(7v) 



dt 
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(196) 



(20) 



(21) 



(22) 



for the spatial equations. Note the reappearance of 7 in the Coriolis and centrifugal forces, 
since 7 is no longer conserved. 

We call attention to the term 



(ft x r) • E 
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in equation (21). It shows that a charged particle may gain or lose inertia even when 
undergoing pure corotation, provided there is a component of E in the direction of fi x r. 
This is obvious when considering particle motion in the nonrotating coordinates, but is easy 
to forget when passing to the rotating frame, where the 'motion' becomes hidden. 

We also see that the Lorentz force takes on a peculiar form in our rotating frame. The 
term E + vxB seems to lend credibility to naming E and B as the electric and magnetic 
fields, but note that we could have also have written it as E + v' x B. 

What seems to be absent in existing literature is the final term in the Lorentz force, 
which can also be written as 

Note that the direction of this force is against (with) the rotation for a particle which is 
gaining (losing) energy, betraying its origin as an inertial 'relativistic mass' effect. 



5. Interpretation 

We have already seen in equation (14) that p and j can be naturally identified with the 
charge and current densities in this frame; we wish to make similar identifications of the 
electric and magnetic fields. The standard technique from special relativity determines the 
electric field by contracting U a with the mixed representation of F: 

E a = (0, E) = F a p U 13 , 

but this doesn't work here due to the nonzero component F t t in equation (20), which is 
also responsible for the (Q x r) • E term in equation (21). (It is interesting to note that it is 
this 'azimuthal' component of E for which E and E always agree.) Clearly, no prescription 
for 'electric' and 'magnetic' fields will reproduce the Lorentz force as we are most familiar 
with it. 
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The root of the problem is seen by considering equation (19b), which says that the 
Lorentz force modifies geodesic motion by 'rotating' the tangent U a to a particle's worldline 
away from its parallel-transported value (along the worldline). These 'rotations,' by (the 
'unprimed' version of) the metric equation (2), are restricted to satisfy the equation 
QjxvU^ dU v jolt = 0. In the flat Minkowski metric of special relativity, these rotations are just 
the directions of Lorentz boosts and spatial rotations. We naturally identify forces which 
produce spatial rotations of U a as magnetic in nature, and those which boost U a as electric 
(for electrically charged particles). In this rotating reference frame, timelike and spacelike 
components are not orthogonal to one another, so boosts and rotations necessarily become 
mixed. Indeed, infinitesimal Lorentz boosts and rotations transform in exactly the same 
way as do E and B, respectively. We conclude from this that any reference frame with 
non-diagonal metric tensor will suffer the loss of precise definitions of electric and/or 
magnetic fields. 

Nowhere is this more apparent than in equation (21), which allows a charged particle's 
inertia to change provided (f2 x r) • E is nonzero, even when the particle is at rest in the 
corotating frame. Conservation of angular momentum demands that this exerts a torque on 
the source of the fields (e.g., the neutron star), which can alter the frequency of rotation. 
This phenomenon could be put to beneficial use in pulsar models. The 'braking index' of 
pulsars, defined as 

n = —. — 

ft 2 

(where the dots denote time derivatives) has, when observable, always been found smaller 
than the canonically predicted value n — 3. (This is predicted as a lower limit when the 
torque is from simple multipole radiation — see, for example, Kaspi et al. 1994 and 
references therein.) Presently, the discrepancy is understood as being due to particle 
outflow along lines of B, or due to more complicated effects. We see, however, that even 
the corotating portion of a magnetosphere, with B • E = 0, can alter the rotation rate by 
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doing work on charged particles (if E • (O x r) 7^ 0). 

It is a pleasure to thank my colleagues Tim Hankins, Jean Eilek, Jim Weatherall, 
David Moffett, and Tracey DeLaney for our weekly pulsar discussions, which motivated the 
present work. Funding for this work was provided by NSF grant AST-9315285. 
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